Quantum tunneling of a single semifluxon in a O-rr Josephson junction 
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We consider a symmetric 0-n Josephson junction of length L, which classically can be in one 
of two degenerate ground states f or J., corresponding to supercurrents circulating clockwise or 
counterclockwise around the 0-tv boundary. When the length L of the junction becomes smaller 
than the Josephson penetration depth Aj, the system can switch from one state to the other due 
to thermal fluctuations or quantum tunneling. We map this problem to the dynamics of a single 
particle in a periodic double well potential and estimate parameters for which macroscopic quantum 
coherence may be observed. We conclude that this system is not very promising to build a qubit 
because (a) it requires very low temperatures to reach the quantum regime, (b) its tiny flux is hard 
to read out and (c) it is very sensitive to the asymmetries between the and n parts of the junction. 
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PACS numbers: 74.50.+r, 75.45.+j, 85.25.Cp 03.65.- 

I. INTRODUCTION 

Josephson junctions with the phase drop of it in the 
ground state (ir JJs)^ are intensively investigated as 
they promise important advantages for Josephson junc- 
tion based electronics 2 -! 3 -, and, in particular, for JJ based 
qubits^*^. Nowadays a variety of technologies allow to 
manufacture such junctions 7 - : 8 ' 9 ' 10 : 11 . 

One can also fabricate so-called 0-ir long Josephson 
junctions (0-7T LJJs ) 12 ' 13 ' 14 ' 15 , i.e., LJJs some parts of 
which behave as junctions and other parts as n junc- 
tions. The most interesting fact about such junctions is 
that a semifluxo n 16 : 17 , i.e. a vortex of supercurrent, car- 
rying one half of the magnetic flux quantum $0 ~ 2.07 x 
10~ 15 Wb, can be formed at the boundaries between the 
and 7r regions provided the JJ is long enough. Classi- 
cally, such a 0-7r LJJ has a degenerate ground state cor- 
responding to cither positive or negative polarity of the 
semifluxon, which we denote as the | or j state, respec- 
tively For these two polarities the circulations of the su- 
percurrents and, therefore, the resulting magnetic fields, 
have different directions. The classical properties of 
scmifluxons are under intense theoretical and experimen- 
tal investigations 1 ^^^^ 2 ^^^^^^^ 2 ^^. 
While the classical properties of scmifluxons (at least for 
systems with few scmifluxons) are more or less under- 
stood, their quantum behavior and their possible appli- 
cations in the quantum domain still have to be studied. 

When the energy barrier separating the two degener- 
ate classical ground states is small enough, the system 
may switch from one state to the other due to thermal 
excitation over the barrier or due to quantum tunnel- 
ing through the barrier. Thermally induced flipping of 
a single semifluxon was already observed 2 ^. Macroscopic 
quantum tunneling (MQT) and macroscopic quantum co- 
herence (MQC) in the system of two coupled scmiflux- 
ons in a 0-7T-0 LJJ was investigated theoreticall y 32 : 33 . In 



Ref. [33| it was also concluded that, unlike a pair of semi- 
fluxons, a single semifluxon in a 0-7T LJJ with moderate 
or large length is always in the classical regime. 



In the present paper, we study a 0-n JJ of length 
L < A j, which admits semifluxon-likc solutions for any 
junction length provided the lengths and critical current 
densities of the and 7r parts are equal. We map the full 
problem to the dynamics of a single particle in a double 
well potential with periodic boundary conditions and es- 
timate the crossover temperature as well as the tunneling 
rate (energy level splitting) between the states t and J.. 
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FIG. 1: (Color online) Schematic drawing of a 0-7T JJ and 
magnetic field profiles in two degenerate ground states. Gray 
curves show the magnetic field profile in an infinite JJ (solid 
line is semifluxon, dashed line is antisemifluxon, $ = ±$o/2). 
Black curves show corresponding solutions in the JJ of finite 
length. In this case the magnetic flux |$| < $o/2. 
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II. MODEL 

We consider a "long" one dimensional Josephson junc- 
tion of length L, see Fig. [T] The Josephson phase n(x, t) 
is a continuous function of the position x (—L/2 < x < 
+L/2) measured along the J J and of time t. We restrict 
ourselves to an undriven dissipationless system. The dy- 
namics of such a system is described by the Lagrangian 
C = K — U , where 



K = Ej 

represents the kinetic energy and 



+ L/2 2 

UJ p — dx 
-L/2 Z 



(1) 



U = E, 



+L/2 



-L/2 



A J^f + 1 ~ cos [m(^: t) + 6{x)] \ dx 



(2) 

is the potential energy. The subscripts x and t denote 
the partial derivatives with respect to position and time, 
respectively. In the above equations the three physical 
parameters are the Josephson energy per unit of junc- 
tion length Ej, the Josephson penetration depth Aj and 
the Josephson plasma frequency tu p . The function 9(x) 
describes the position of and 7r regions along the junc- 
tion. It is zero along regions and is equal to tt along tt 
regions. 

Applying the Euler-Lagrange formalism to our La- 
grangian, we find that on the classical level the dynam- 
ics of the Josephson phase is described by the time- 
dependent sine-Gordon equation^ 

\ 2 j^ xx - w" 2 /i tt - sin[/i(a;, t) + 6(x)\ = 0. (3) 



III. A FRACTIONAL VORTEX IN A 
SYMMETRIC 0-tt JUNCTION 

Let us consider a symmetric 0-7r LJJ of the length L 
(—L/2 < x < +L/2) with a 0-7r boundary at x = 0, see 
Fig. [T] In this case 9(x) is a step function 



9(x) 



0, -L/2 < x < 0; 

7T, < x < +L/2. 



(4) 



Classically, for any L the trivial stationary solutions 
fi(x) = or ^(x) = 7r are unstable and correspond to the 
Josephson energy maximum C/ max = EjL, sec Eq. ©. 
The stable stationary solutions depend on the length L 
of the system. 

If L ^> A j, the ground state of this system is a single 
scmifluxoni 12 i 16 ' 17 i 22 Such a semifluxon may have posi- 
tive or negative polarity that corresponds to two classi- 
cal degenerated states | and | 16 i 22 , see the gray curves 
in Fig. [T] In this limit the energy barrier separating two 
classical states is very large and the system is always in 
the classical regime^. Therefore, in the present paper 
we investigate 0-n LJJs with L < Aj. 



For L < \j the stationary stable solutions of Eq. ([3]) 
can be approximated b y 34 i 35 



7T 

± 2 ± 



<L-\x\) 
2A 2 



(5) 



and correspond to the energy minima. Magnetic field 
profiles are shown in Fig. [T] by black lines. Substituting 
the stable and the unstable stationary solutions for [i into 
Eq. ([2]), we find that the barrier between two stationary 
stable solutions is given by 



U (£) w Ej\j 



£3 
24' 



(6) 



where I = L/Xj is the normalized J J length. In each of 
these ground states the flux in the junction is 

$0 -A M = p. [fx(L/2) fx(-L/2)} = ±* £. (7) 



<I> 



2tt 



2tt 



'8tt 



The main assumption of this paper is that for suffi- 
ciently short LJJ and sufficiently low energies the time 
dependent phase n(x, t) can be approximated by (see Ap- 
pendix [AJ 



^{x,t)=Q{t) + \sHQ{t)f 



A 2 , 



(8) 



i.e. it moves along a collective coordinate Q. Note that 
the ansatz ([5]) satisfies zero magnetic field boundary con- 
ditions fi x (±L/2) = 0. The values of Q = ±ir/2 cor- 
respond to the two distinct classical ground states ((5]), 
whereas the values Q = or Q = ir correspond to unsta- 
ble solutions /i = and fi = -k. 

By substituting the ansatz ((8]) into Eq. ([2]) we calculate 
the potential energy of the system as a function of Q and 
obtain 



C/(Q)= J B J A,^(l--sin 2 Q 



(9) 



Clearly, the barrier between the two classical ground 
states is given by Eq. ^ which simultaneously is the 
amplitude of the cos-like potential. 

In essence, our ansatz ([5]) maps the full problem to 
the motion of a fictitious particle in the one-dimensional 
potential U(Q) ©. 

To find the mass of this particle we substitute the 
ansatz ([5]) into Eq. (TTJ) and find 



K(Q) 



Ejt Q 2 



(10) 



Since K(Q) = MQ 2 /2, the mass of our fictitious particle 
is given by 



M 



E 7 X 



[kg • m 2 



(11) 



Note, that actually M is not a mass, but a moment of 
inertia, since our coordinate Q is not a position, but a 
dimcnsionlcss phase. 
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A. Estimate of quantum-to-classical crossover 

For a harmonic oscillator with mass M and frequency 
luq the (square of the) width of the probability distribu- 
tion in the ground state is determined by 



(AQ 2 ) = 



Mujq 



(12) 



Our potential U(Q), Eq. is not parabolic, but wc 
may approximate each potential well with a minimum at 
Qo = ±7r/2 by a parabolic potential, i.e., 



U(Q) « EjXj 



24 



24 



(Q - Qo) 2 



(13) 



which together with the mass M given by Eq. (fTTj) de- 
scribes a harmonic oscillator with the "spring constant" 
k = EjXj£ 3 /12 and the eigenfrcquency 



LUQ 



(14) 



Using our expressions for the mass M (fill) , the frc 
quency loq (|14p . and the harmonic oscillator width JT - 
we obtain 



(AQ 2 ) = 



EjXj £ 2 



(15) 



as an estimate for the spread of a wave function in each of 
the wells of the potential U (Q). Quantum effects are no- 
ticeable when the wave function in the left well overlaps 
with the wave function in the right well. This overlap 
should be appreciable, but not too large since otherwise 
the two states will not be distinguishable. For a rough es- 
timate we take (AQ 2 ) > 0.1Q 2 , as criterion for quantum 
behavior. Thus quantum effects will dominate if 

(AQ 2 ) _ 8^3 fruj p 1 



ir 2 EjXjP 
32^ h I [i d ! 1 



1 % V Cw 2 P 
where we took into account the definitions 



£0.1, 



Xj = 



$0 



2-K(j, d'j c 



E.J = 



2tt 



(16) 



(17) 



In Eq. (|17[) . fiod' is the inductance per square of the su- 
perconducting electrodes (fia is the permeability of vac- 
uum, d! « 2Ai, Xl is the London penetration depth), j c 
is the critical current density of the LJJ, C is the capac- 
itance of the LJJ per unit of area, and w is the LJJ's 
width. 

For typical parameters Xl = 90 nm, w = 1 fiia and C = 
4.1/iiF/cm 2 (Hypres^ technology with j c = 100 A/cm 2 ) 
Eq. ((T6|) reads 



(AQ 2 
Qo 



3.2 x 10- 3 r 2 > 0.1 



(18) 



Thus, quantum effects start to play a role for £ < 0.18. 
Note, that according to Eq. (flB|) the classical-to-quantum 
crossover length does not depend on j c . 

Using definitions (|17p in terms of physical parameters 
of the LJJ, we can express the inertial mass M (|lip as 



M 



<jC 



\/fJ-od'jc 



2n 



1.3 x 10" 4 £m e A 2 ,. 



(19) 



This means that a single electron moving around the 
whole JJ (circulating around 0-ir boundary) has much 
larger moment of inertia than our semifluxon. 

To estimate the classical-to-quantum crossover tem- 
perature, we compare the thermal escape rate and the 
scmiclassical expression for the quantum mechanical en- 
ergy level splitting. The thermal escape rate is cx 
exp(—Uo/kBT), while in the semiclassical limit the en- 
ergy level splitting is cx exp(— 4?7 /^o), see Eq. (|B9[) . 
Neglecting prefactors we conclude that at the tempera- 
ture 



4fc^ 



fku p £ 
8V3k B 



(20) 



quantum tunneling and thermal escape have the same 
rate, while at lower temperature quantum tunneling 
dominates. For the parameters chosen above and £ — 
0.18 this gives T x w 27 mK. 

Another specific temperature T* is defined by compar- 
ing k B T with the barrier height. For T > T* the states j 
and { cannot be distinguished, i.e., even in the classical 
regime read out cannot be properly organized. For our 
Joscphson junction the temperature T* is given by 



EjXj £ 3 



24 



(21) 



For £ = 0.18 we obtain T* w 218 mK. Note that T* can 
become smaller than T x for small values of I. For such 
parameters, however, we have XJ§ < hw and we are not 
in the semiclassical limit used to derive Eq. (|20l) . 



B. Energy level splitting 

The stationary Schrodinger equation for our collective 
coordinate Q reads 



if 



2M dQ 2 



U{Q) 



xjj = EtJj. 



(22) 



Since in our configuration the Josephson phases [i + 2ix 
and fi cannot be distinguished, fi (and therefore also Q) 
is only defined modulo 2tt. Therefore, we supplement 
Eq. ([22|) with periodic boundary conditions ip(Q + 2ir) = 
V(Q). 

The main purpose of our quantum mechanical calcula- 
tions is to investigate quantum tunneling between the two 
classical ground states. More precisely, we want to ap- 
proximate our Joscphson junction by a quantum mechan- 
ical two-level system which may be used to implement a 
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qubit. In such a two-level system the dynamics of the 
two lowest levels is decoupled from the other levels. This 
requires that the two lowest energy eigenvalues E\ and 
E-2 of the Schrodinger equation (f2"2")) are well-separated 
from the other energy eigenvalues. 

Furthermore, we want to observe coherent oscillations 
between the left and the right well of the potential. This 
is only possible if there are superpositions of the ground 
state -01 (Q) an( l the nrs t excited state ip2(Q) of the 
Schrodinger equation (f2"2")l which are sufficiently well lo- 
calized in the left and in the right potential well. 

We expect that we can approximate our system by a 
quantum mechanical two-level system, if the two classi- 
cal (degenerate) ground states are separated from each 
other by a sufficiently high energy barrier. In our sys- 
tem, the barrier height depends on the scaled length £ of 
the Joscphson junction, see Eq. (|6|). To gain more insight, 
we have solved our Schrodinger equation (|22p numerically 
for different values of £. The parameters of the Josephson 
junction are given in Sec. IIII Al above Eq. I|18p. 
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FIG. 2: The two lowest energy levels Ei and E 2 , AE = E 2 - 
Ei and the barrier height Uq (all given in mK) as a function 
of I for typical parameters given in the text. 

Fig. [U shows the lowest energy levels E\ and E 2 , the 
energy level splitting AE = E 2 — E\, and the barrier 
height Uq of the potential as a function of the normal- 
ized JJ length £, When we reduce £, the lowest energy 
eigenvalues split. For I < 0.1 the energy level E 2 is lifted 
above Uo so that we have only one energy level below Uo. 

Fig.[3]shows the eigenfunctions ipi(Q) and tp 2 (Q) corre- 
sponding to the eigenvalues E\ and E 2 . We have chosen 
two values for i which lead to a qualitatively different 
behavior. 

For I = 0.13 (Fig. [3^,) we have a quantum mechanical 
two-level system as discussed above: The two lowest en- 
ergy eigenvalues E\ and E 2 are well-separated from the 
others and ["0i(Q) ± ip 2 (Q)] /v2 is well-localized in one 
of the potential wells. 

For I = 0.09 (Fig. [3Jd) our simple picture of a quan- 
tum mechanical two-level system fails: The difference be- 
tween Ei and E 2 is large and only one energy eigenvalue 
is below Uq. Furthermore, the two lowest eigenvalues 
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FIG. 3: Solutions of the Schrodinger equation (|22|l for (a) 
£ = 0.13, (b) ^ = 0.09. 



are not well-separated from the higher energy eigenval- 
ues (not shown in Fig. [3Jd). The ground state ijJi(Q) is 
completely delocalized and \ip\{Q)\ 2 does not have pro- 
nounced maxima at the minima of the potential. Obvi- 
ously, for £ = 0.09 we cannot approximate our system by 
a two-level system. 

As expected, only for sufficiently large values of £, i.e., 
sufficiently high energy barriers, our Josephson junction 
behaves like a two-level system. However, we have to 
take into account that thermal excitations will destroy 
quantum coherence if the energy level splitting AE be- 
comes too small. Present technology restricts AE to val- 
ues larger than 20 ... 30 mK — the temperature achiev- 
able in modern dilution refrigerators. According to Fig. [5] 
this requires £ < 0.09. For such values of £, we are not in 
the limit of a two-level system, see Fig. [3b. 



IV. DISCUSSION 

The results for the parameters used in Fig.[5]and Fig. [3] 
are not very promising: cither we are not in the limit of 
a two-level system or we need temperatures significantly 
lower than 20 ... 30 mK to avoid thermal excitations. We 
can try to improve the situation by using different param- 
eters for the Josephson junction. 

In Appendix IB II we show that the Schrodinger equa- 
tion (|2"2"|) is equivalent to the Mathicu equation (|B4[) with 
only one scaled system parameter h and the scaled energy 
e, see Eqs. (|B2| and ([B3)) . Using the definitions Ifl?]). we 
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(a) single semifluxon 



I 


T* 


T x 


AEf 2 


AEi2 m 


AE^ m 


0.18 


218 


27.0 


0.3 


0.2 


89.6 


0.13 


83 


19.5 


7.2 


5.3 


40.8 


0.10 


38 


15.0 


24.4 


15.5 


18.9 


0.09 


28 


13.5 


32.0 


19.9 


18.8 



(b) AFM semifluxon molecule 



5a 


T* 


T x 


AE? 2 


AE?% m 


AEgp 


0.020 


566 


58 


0.2 


0.2 


265 


0.015 


318 


51 


3.4 


2.8 


183 


0.010 


141 


41 


36.5 


22.9 


124 



TABLE I: The crossover temperatures T* and T x , the energy 
level splitting AE\2 = E2 — E\, and the energy difference 
AB23 = Es — E2 between the first and the second excited 
state for a single semifluxon (a) and a semifluxon molecule (b) 
for different values of the scaled lengths £ and Sa. All values 
are given in mK. The numbers for the semiclassical energy 
level splitting Ai?i2 are based on Eq. {B9) (a) and Eq. (fB10)) 
(b), whereas AS"™ and A-EJj™ were calculated numerically. 
For A_Bi2 <C AE23 the systems are good two-level systems. 



find that the parameter h is proportional to (j c w) 2 and 
that the scaling factor between the energy E and the 
scaled energy e is proportional to w^ 1 and does not de- 
pend on j c . If we reduce the width w by a factor of 
10 and simultaneously increase j c by a factor of 10, we 
can increase the energy scale by a factor of 10 without 
changing the generic behavior of the system (h does not 
change). For the parameters used in Fig. we could 
increase the energy level splitting from 5.3 mK to 53 mK 
and would still be in the limit of a two-level system. 

To read out the state of a semifluxon one can use a 
SQUID situated just in front of the 0-ir boundary. Note, 
however, that the two states carry the flux <f> <C &o/2. 
For small £ the flux is given by Eq. ([7]) and is equal to 
~ 10~ 3 $ for £ = 0.18. Needless to say that it is rather 
difficult to measure such a small flux accurately. 

Another problem is a high sensitivity to the parameter 
spread. All calculations done so far are for a symmetric 
0-7T junction, i.e., for the ideal case where the critical 
current densities and lengths of the and it parts are 
equal. In experiments, though, a small asymmetry is 
always present. For ^ < 1 this may easily lead to the 
situation where the classical ground state is just fj, = 
(when the part is longer or has a higher j c ) or fj, = n 
(when the 7r part is longer or has a higher j c ) which is 
not doubly degenerate. For £ <C 1 the length asymmetry 
should not exceed^ £ 3 /48, which is - 1.2 x 10" 4 for £ = 
0.18. Such precise control of LJJ length is not feasible. 

Finally we want to compare our results to the results 
we obtained for an AFM semifluxon molecule in a 0- 
7T-0 junction 3 - 3 .. In such a system the main parameter 
which controls the height of the energy barrier (similar 
to £ in this paper) is the length a of the ir region. It 
turns out that to approach the quantum regime a should 



be just a bit above a crossover length a c = ttXj/2, i.e., 
a = (rr/2 + Sa)Xj. The corresponding key numbers are 
given in Tab.[T)D. According to the Tab.Uthe semifluxon 
molecule in a 0-7T-0 junction is more suitable to build 
qubits. The crossover temperatures for Sa = 0.01 can be 
reached with modern dilution refrigerators. 

V. CONCLUSIONS 

We have presented a simple quantum theory of a short 
0-7T JJ. By introducing a collective coordinate, the full 
problem was reduced to a dynamics of a single particle 
in a one-dimensional potential. 

It was found that the quantum regime is reached when 
the normalized length £ of the junction is smaller than 
0.18. The classical-to-quantum crossover temperature T x 
is about 20 mK for the case when we still have a good two- 
level system and the semiclassical approximation (used to 
define T x ) is still valid. This is on the limit of modern 
dilutions refrigerators. Moreover, even for I < 0.18 the 
flux, which should be read out to determine the state 
of the system, is < 10~ 3< E>o, which is rather difficult to 
detect even using SQUIDs. The system is also very sen- 
sitive to asymmetries in the and ir parts, which should 
not exceed 10 ~ 4 . 

It is quite interesting that another system, the so-called 
"d-dot"—, is equivalent to the 0-7T LJJ of the present 
paper. A "d-dot" is the squared island of a d-wave su- 
perconductor embedded into an s-wave subspace. Es- 
sentially it is an annular 0-7T-0-7T LJJ with the length of 
each region equal to a. The two ground states arc tltl 
and I tit- in iac t, the four regions of a d-dot were used 
just because of the topological limitations. One can use 
symmetry arguments to conclude that such an annular 
0-7r-0-7r JJ is equivalent in terms of ground states and 
cigenmodes to an annular 0-7T JJ with only two facets of 
length a, which, in turn, is equivalent to a linear (open 
ended) 0-7r J J with the total length a. In contrast to 
Ref. 1371 . we conclude that such a system is not a very 
promising candidate to build qubits. Better figures are 
shown by a two-scmifluxon molecule in a 0-7T-0 JJ with 
tunneling between "fj and || state a 32 i 33 . 
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APPENDIX A: DERIVATION OF COLLECTIVE 
COORDINATE 

In this appendix we motivate the use of the collective 
coordinate Q(t), that is, the assumption that for suffi- 
ciently short LJJs and sufficiently low energies we can 
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approximate the Josephson phase t) by 

x (L — \x\ 



H{x,t) ss Q{t)+sm[Q(t)}- 



2A5 



(Al) 



where Q(t) describes the motion of a single particle of 
mass 



M 



Ej\j 



in the potential 



U(Q)=EjXj£[l-—sm' 



(A2) 



(A3) 



In order to keep the notation simple, we use the scaled 
length £ = L/Xj of the Josephson junction, the scaled 
time r = ujpt and the scaled position f = x/Xj in the 
rest of this appendix. 



1. Fourier expansion 

We start from the Lagrangian 

+t ' 2 / n 2 \ 2 

C = EjXj >^ 9 ^ ^ dfi 



2\dr 



2\d£ 



- [1 - cos (/x + 0)] \ df (A4) 



-1/2 



for the Josephson phase /i(f,r) and expand /x(f,-r) into 
a Fourier series which takes into account the boundary 
conditions /j,^(—£/2,t) = ^(£/2,t) = 0. The most gen- 
eral Fourier series that takes into account these boundary 
conditions is 



t (f,r) = Q(r) + £ ? W(r) cos ^f) 

OO 



ra=0 



/2n+l . 
sm I — - — 7rf 



(A5) 



In the following calculations we make various approx- 
imations for small I. A detailed analysis shows that the 
cosine terms in this Fourier scries only lead to higher- 
order corrections in £ 2 which are to a large extent not 
consistent with the approximations we are going to make. 

In order to keep the calculations brief, we therefore 
substitute the Fourier series 



H&t) = Q(T) + V2j2ln(T)sm 



2n + l 



£ 



< (A6) 



into the Lagrangian (|A4[) . The factor \[2 in front of the 
sum was chosen for convenience, see the kinetic energy 
in Eq. (|A11[) . For short LJJs. non- uniform excitations 



require large energies since for small values of £ and sig- 
nificantly large non-uniform excitations the term dfi/dt; 
in the potential energy of the Lagrangian (|A4[) becomes 
large. Therefore, we expect that the Fourier coefficients 
q n {T~) are small for low energies. Since we are interested 
in the low-energy behavior of the system, we only take 
into account terms in the Lagrangian that are quadratic 
in q n (r) and and neglect highcr-oder terms. For the 
uniform part Q(t), however, we make no assumption. 
A similar approach was used by Koyama et alM- for a 
square-shaped closed 0-7r junction (d-dot). 

Due to the orthogonality of the sine functions in the 
Fourier series (|A6|) the only non-trivial integral we have 
to calculate is 



+t/2 



cos [M£,t) + 0(f)] 



(A7) 



Here we make use of our assumption that the non- 
uniform part 5(f, t) = t) — Q(t) of the Josephson 
phase is small and rewrite cos [/Lt(f , t) + 0(f)] in the form 



cos W,t) + 0(f)] 



cos Q(t) cos J(f , t) cos 0(f) 
— sin Q(t) sin 5(f , r) cos 0(f) 
cos Q(t) [1-<5 2 (£,t)/2] cos 0(f) 
-sinQ(r) 5((,t) cos 0(f), (A8) 



where we have used sin 0(f) = 0. 

The integral over cos 0(f) alone is zero since the 0- 
region and tt region have the same length. We are there- 
fore left with the two integrals 



6(£,t) cos 0(f) df 



(A9) 



-1/2 



and 



5 2 (f,r) cos 0(f) df. 



(A10) 



Due to the different sign of cos 0(f) in the 0- and in the 
7r region we cannot use the orthogonality relations of the 
sine functions in Eq. (|A6[) . Evaluating the integrals is 
nevertheless straightforward as the second integral van- 
ishes for symmetry reasons. 

Within our approximation we finally arrive at the ap- 
proximate Lagrangian 



L « E.jX 



11 i 

?&+~51fi-V{Q,{q n }) , (AH) 



n=0 



where the dot denotes the derivative with respect to r 



7 



and where the potential energy V (Q, {q n }) is given by 

2 00 



n=0 



.^ sinQ y (A12) 



2. Harmonic oscillators for q n 

For fixed values of Q the potential V (Q, {q n }) is har- 
monic with respect to {q n }- In order to get more insight, 
we rewrite it in the form 



V (Q, {qn}) = V (Q) + ^ >T(2n + l) 2 [q n - q n {Q)f . 

n=o 

(A13) 

By comparing Eqs. (|A12|) and (|Al3j) we find that we 
have to choose 

2^/2/' 2 

g"(Q)= s/o I iv« sin Q (A14) 



7r 3 (2n + l) 3 



and Vb(Q) is given by 



Vo(Q) = l-^f^(2n+lff n (Q) 



n=0 



1 - ^ sin Q, 



where in the last step we have used^ 

1 

(2n + l) 4 



00 4 

El TT* 
(2n + l) 4 = 96' 



(A15) 



(A16) 



For fixed values of Q this potential describes a sum of 
independent harmonic oscillators with (scaled) frequen- 
cies 

fl n = (2n | 1)7r , 71 = 0,1,2,3.... (A17) 



3. Stationary solution 

The stationary solutions of the full potential 
V(Q, {q n }) follow from 

71— 

(Al8a) 

^ = ^(2n+lf[q n ~q n (Q)]=0. (A18b) 

From these conditions we immediately find 

sinQocosQo = 0, g£ = q n (Qo)- (A19) 



We therefore have two classes of stationary solutions: 
sinQ o = 0, Qa = kTT (k = 0, ±1, ±2, . . . ), 

<7°=fe(Qo) = (A20) 

and 

cosQ = 0, Q = Ti"/ 2 + k7r (& = 0,±1,±2, . . .), 
(-l) fc 2V2^ 2 



Qn =q n {Qa) = 



7r 3 (2n + l) 3 



(A21) 



A stability analysis shows that the stationary solutions 
corresponding to Eq. (|A20[) are unstable whereas the the 
stationary solutions corresponding to Eq. (|A21j) are sta- 
ble. 

The first class of stationary solutions, Eq. (|A20|) , ob- 
viously corresponds to the uniform stationary Joseph- 
son phases fJ^(0 = kir. When we substitute Eq. (|A21[) 
into the into the Fourier series (|A6|) we obtain the non- 
uniform stationary Joscphson phases 

o,m * , (-l) fc 4^ 2 

Z 7T° 



n=0 



(2n+l) 



(2n + l e 
sin I — - — ire, 



| + tor + ^£('-|£l)> (A22) 



where we have used^ 

E 3 sin [(2n + l)x] tt 
n (2n+l) 3 8 

n— v 1 



x(jT — \x\), —IX < X < TT. 



(A23) 



4. Reduction to a one-dimensional system 



a. Harmonic oscillators at low energies 



In Sec. IA 21 we found that for fixed values of Q the po- 
tential V(Q, {q n }) describes a sum of independent har- 
monic oscillators with frequencies Q n = (2n+ 1)tt/£,ti = 
0, 1, 2, 3, ... . For decreasing £ the harmonic potential be- 
comes steeper and steeper whereas the barrier height in 
the potential Vq(Q) decreases, see Eq. (|A15[) . Therefore, 
we expect that for small £ and low energies the dynam- 
ics of the full system is restricted to small fluctuations 
around q n (Q) and we may approximate <7n(r) by 



q n (r) ~ q n (Q(r)) ■ 



(A24) 



When we substitute this approximation into the La- 
grangian (|A11|) we obtain 



A \q 2 



dq n 
dQ 



V (Q)\. (A25) 
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The terms with the derivatives of q n (Q) are of the order 
£ A whereas in Vq{Q) only the terms ~ £ 2 are taken into 
account. Therefore, we decided to neglect these terms 
and finally arrive at the effective Lagrangian 



£ c ff f» EjXjt 



Q 2 - V {Q) 



which describes the motion of a particle of mass 
M = ^£ 

in the potential 

U(Q) = EjXjl (l - ^ 



sin Q 



(A26) 



(A27) 



(A28) 



The dynamics of a particle in this potential can be char- 
acterized by the barrier height 



of the potential U(Q) and the frequency 



w = 



2\/3 



(A29) 



(A30) 



for small oscillations around the minima of the potential 
U(Q). 

The Fourier series (|A6|) connects Q{t) and q n {T~) 
with the Joscphson phase /z(£,r). When we substitute 
Eq. (|A14[) into the Fourier series (|A6|) the corresponding 
approximation for the Joscphson phase r) reads 



U 2 

<3( T ) + sin g(r) 



E 



i 



^ n (2n + I) 3 

n— 7 



Sill 



2n+ 1 



Q(r) + - sin Q(r) • t(£-\£\) 



(A31) 



where in the last step we have used Eq. (|A23[) to evaluate 
the sum. Equation (|A31[) is the central result of this 
appendix and is used as the starting point in Sec. lIIIi sec 
Eq. ©. 

Finally we want to mention that our results are in 
agreement with the assumption of Sec. I A 1 1 that the non- 
uniform part of the Joscphson phase is small. 



b. Quantum mechanical considerations 

So far our discussions were purely classical. In the 
main part of the paper, however, we use Q(t) as a col- 
lective coordinate to investigate quantum tunneling in 
Joscphson junctions. We now briefly estimate under 



which conditions our collective coordinate Q(t) is use- 
ful in the quantum limit of the system. 

We have several energy scales in the system: the ther- 
mal energy ksT, the barrier height Uo, Eq. (|A29|) . of the 
potential U(Q), and the excitations hw p D. n of the har- 
monic oscillators discussed in the Sec. IA2I Using Q(t) 
as a collective coordinate in the quantum regime obvi- 
ously only makes sense if these harmonic oscillators are 
not excited. Therefore ksT and Uo have to be much 
smaller than hw p ^l n . Since we are mainly interested in 
quantum tunneling in the potential U(Q), ksT has to 
be smaller than Uo, see also the discussion in Sec. MI Al 
The condition that the harmonic oscillators are not ex- 
cited therefore reduces to 



U -C hujpilo 



or 



4 2^hw p 
EjXj ■ 



(A32) 



(A33) 



Furthermore, we have to restrict ourselves to £ 2 <C 1 since 
various approximations in this appendix are expansions 
in £ 2 . If these condition are fulfilled we can assume that 
all harmonic oscillators are in the ground state. 



c. Semiclassical limit 

Our first attempt to identify a collective coordinate 
was based on the intuitive picture of harmonic oscillators 
which are not excited for small values of £. In a second 
independent approach we determine a trajectory which 
can be used to describe quantum tunneling in the full 
system in the semiclassical limit. 

A standard method to derive semiclassical expressions 
for quantum tunneling are instanton a 39 ' 40 i 41 ' 42 . In this 
method Feynman path integrals are extended to imag- 
inary times. For simple systems where the Lagrangian 
can be written as the difference between kinetic energy 
and potential energy this is equivalent to path integrals in 
real time for the inverted potential. As in the ordinary 
Feynman path integral formalism in quantum mechan- 
ics, semiclassical quantum mechanics is based on classical 
trajectories and small fluctuations around these trajecto- 
ries. 

We use the instanton technique only to determine the 
classical trajectory that governs quantum tunneling in 
the semiclassical limit. As in our first approach this tra- 
jectory is then used to define a collective coordinate. In 
this appendix we do not use instantons to derive semiclas- 
sical expressions for the full system. This would require 
to take into account fluctuations around this trajectory. 

In our simplified approach we have to find the classical 
trajectory that starts with zero kinetic energy at one of 
the maxima of the inverted potential and ends with zero 
kinetic energy at the other maximum. Due to the sym- 
metry of the potential V(Q, {q n }) it is sufficient to find 
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the trajectory which starts at 

Q = 0, q n = (A34) 
(minimum of the inverted potential) at r = and reaches 



<l« 



7r 3 (2n+ l) a 



(A35) 



(maximum of the inverted potential) for r — > oo. 

If we assume that we can express q n as a function of Q, 
that is, g n (r) = q n (Q(i~)), we can rewrite the equations 
of motion of a particle moving in the inverted potential 
— V(Q, {q n }) in the form 



W dQ 2 dQ dQ 
Energy conservation requires 



ay 



0. 



(A36a) 
(A36b) 



n=0 



fdq n \ 2 
\dQJ 



-V(Q,{q n }) 



24 



(A37) 



which can be used to express Q 2 as a function of Q. Sub- 
stituting this relation into Eq. (|A36b[) we obtain, together 
with the boundary conditions Eqs. (|A34[) and (|A35[) , a 
system of coupled nonlinear differential equations for the 
trajectory we are looking for. 

Instead of solving the differential equations for q n (Q) 
approximately we use the results we obtained so far. 
These results suggest that q n (Q) is approximately given 
by qn{Q) ~ <7n(Q)- We verify this by substituting 



q n (Q) = q n {Q) + i 6 Vu{Q) 



(A38) 



into Eqs. (|A36bl) and ||Al7|). W e assume that ij n (Q) is of 
the order one. From Eq. (|A37[) we obtain 



ri2 



12 



cos 2 Q + 0(£ 6 ). 



Furthermore, we use 



BV £ 2 

— = --sinQcosQ + 0(0, 
F)V 

= ir 2 £ 4 (2n+l) 2 r, n (Q). 



(A39) 

(A40a) 
(A40b) 



If we take into account that q n {Q) is of the order i 2 [see 
Eq. (IA14[) ] the lowest-order terms of Eq. (|A36b[) read 

£ 2 d dn 

-cosQ-^cosQ^-7r 2 ^ 4 (2n + l)V(QH0 (A41) 
12 dQ dQ 

and we immediately find 
V2 



Vn(Q) 



37r 5 (2n+ l) 5 



sin Q cos Q. 



(A42) 



This result confirms that ij n (Q) is indeed of the order 
one. We can therefore approximate the trajectory along 
which we have to calculate the Euclidian action So by 



q n {Q) « q n (Q)- 



(A43) 



With the help of Eq. (|A39|) and after omitting a con- 
stant term in the potential energy we can calculate the 
Euclidian action along this trajectory without knowing 
Q(t) as follows: 



So 

T 



EjXjl 



tt/2 

EjXji f 1 



24 



cos 2 Q 



d.T 



-tt/2 
tt/2 



Q 



EjXjt 



1 EjXj 2 



1 • £ 2 
^Q 2 + 24 cos 2 Q 



cosQ dQ 



dQ 



V3 ^ 



-tt/2 
l 2 = 



12 

4^o 

tlWo ' 



(A44) 



Here we only have taken into account terms of the order 
£ 2 . In the last step we have used the definitions of Uo 
and uj , Eqs. (|A"29|) and (|A"30|> . 

It is important to note, that Eq. (|A44[) can be read in 
two ways: According to our calculations, it describes the 
approximate Euclidian action along the trajectory de- 
fined by q n (T) = q n {Q{T~)). Obviously, it also describes 
the exact Euclidian action of a particle moving in the 
potential U(Q) = EjXji 3 cos 2 Q/24. Therefore, in the 
scmiclassical limit we obtain the correct exponent for the 
energy splitting when we use Q(t) as a collective coordi- 
nate to approximate the full system by a system with a 
single degree of freedom. 

In the semiclassical limit the exponent So/h becomes 
large. Taking into account that various approximations 
in this appendix are expansions in £ 2 , the scmiclassical 
limit is characterized by 



Hujp 

eJx. 



«.£ 2 < 1 



(A45) 



within the framework of this appendix. If the parame- 
ters of the Josephson junction fulfill this condition, our 
collective coordinate Q(t) seems to be a useful concept 
to describe quantum tunneling in short LJJs. 



5. Summary 

In this appendix we have shown that for appropri- 
ate parameters, the Josephson phase fi(x,t) can approx- 
imately be described by a collective coordinate Q(t). We 
have expanded the Josephson phase fi(x, t) into a Fourier 
series and used two approaches to reduce the many de- 
grees of freedom of the system to one degree of freedom. 
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In the first approach we used the intuitive picture of 
harmonic oscillators which are not excited for small t. In 
the second approach we determined the trajectory which 
allows us to calculate the exponent of the semiclassical 
expression for the energy splitting. In lowest order in C 2 
we obtained the same collective coordinate Q(t) for both 
approaches. 



APPENDIX B: ANALYTICAL RESULTS 



As discussed in Sec. lIIII and Appendix [Al our collective 
coordinate Q describes the dynamics of a single particle 
of mass M [see Eq. (|A27|) ] in a onc-dimcnsional poten- 
tial U(Q) [see Eq. (|A28j) ]. The corresponding stationary 
Schrodinger equation reads 



2M dQ 2 



U{Q) 



HQ) = Ei/)(Q). (Bl) 



Since we have U (Q) oc cos 2 Q this is essentially the well- 
known Mathicu equation^. 



1. Mathieu equation 

When we take into account the expressions for M, 
Eq. P57|) . and U(Q), Eq. (|M5]) . and introduce the di- 
mensionless parameter 



EjXj 



2 OA 



hujp J 48 \ huo J 
and the scaled energy 



e = 2£ 



EjXjY E 96h E 



hw p ) EjXj e 3 Ej\j 

we obtain the standard form 

ip"(Q) + [e-2h cos(2g)] V(Q) = 



(B2) 



(B3) 



(B4) 



of the Mathieu equation^ 3 -. Here we have omitted an 
unimportant ^-dependent term in the potential energy 

Since the Josephson phases /i(x, t) + 27r and (i(x,t) 
cannot be distinguished, fi(x, t) (and therefore also Q(t)) 
is only defined modulo 2n. Therefore, wc supplement 
Eqs. (|B1[) and (|B4[1 with periodic boundary conditions 
yj(Q + 2vr) = tp(Q). 

From the definition of the parameters h [see Eq. (|B2[) ] 
and e [see Eq. (|B3[) ]. we find that the quantum mechani- 
cal behavior of the system is governed by three parame- 
ters: The purely geometrical scaled length £, the energy 
scale EjXj of the system, and the dimensionless parame- 
ter hw p /(EjXj). We want to emphasize that in Eq. (|B4p 
the only parameter is h, since e is the eigenvalue that we 
want to determine. Furthermore, h does not depend on 
the energy scale EjXj of the system. 



For a given value of h, Eq. (|B4|) has it periodic and 
2-7r-periodic solutions only for characteristic values of e. 
These solutions are either even or odd. Other solutions 
of Eq. dH do not fulfill %j)(Q + 2tt) = HQ). Following 
the notation of Abramowitz and Stegun^, we denote the 
even solutions by ce„(Q, h) and their characteristic values 
by a n (h) (n = 0, 1, 2, 3, . . . ). Similarly, the odd solutions 
and their characteristic values are denoted by sc n (Q,h) 
and b n (h) (n = 1,2,3,... ), respectively. A closer anal- 
ysis shows that the scaled ground-state energy is given 
by e = do whereas the scaled energy of the first excited 
state is given by £ = b\. The corresponding (unnor- 
malized) eigenfunctions are ceo(Q,h) (symmetric and ir 
periodic) and sei(Q, h) (antisymmetric and 27r-periodic), 
respectively. 



2. Approximate results 

Approximations for the characteristic values a n and b n 
are available for h <C 1 and /i > 1, see Ref . We can 
use these approximations to find approximations for the 
difference Ae = £2 — £1 = b\ — ao between the two lowest 
eigenvalues of Eq. (|B4[) . 



a. Low energy barrier (h <C 1) 

For (low energy barrier) we find 

3. 



As 



l-h + -h* 



(B5) 



Using the definitions of h and e, Eqs. (|B2[) and (|B3[) . 
we may therefore approximate the difference AE be- 
tween the two lowest eigenvalues of the Schrodinger equa- 
tion flBl]) by 



EjXj ( hu. 



AE 



2t 



EjX„ 



if £ satisfies 



e 4 < 48 



EjX 





'EjXj 


48 I 




'EjXj 












t) 2 





(B6) 



(B7) 



The lowest-order term corresponds to h = 0, that is, no 
energy barrier. 



6. High energy barrier (h 3> 1) 

For ft > 1 (high energy barrier) we find the exponen- 
tially small difference 



,3/4 -AVh 



(B8) 
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Substituting the expressions for h [Eq. (|B2|) ] . Uq 
[Eq. (|A29]) ]. and uj [Eq. (pOO]) ]. we obtain from Eq. (jB8|) 



A£a;8liwo 



cxp 



_t/o_ 



(B9) 



as the the semiclassical expression for the energy level 
splitting for a particle in a periodic cos(2Q) potential. 
Note, that the exponent is the Euclidian action S$/h [see 
Eq. (|A44j) ] that we calculated in Sec. IA4cl 

When we compare this result to the well-known ex- 
pression 



AE. 



quartic 



, 2U 
>hw \l — — exp 
irnuo 



16 _C7o_ 

3 Hluq 



(BIO) 



for a particle in a quartic double-well potential (see for 
example Ref. |4jjj_2j) , wc immediately see, that these two 
expressions only differ by a numerical factor of the order 
one (4 vs. 16/3) in the exponent. This confirms that the 
asymptotic expression (|B9[) is indeed the semiclassical 
limit for the energy level splitting. Note, that due to 
our periodic boundary conditions, "tunneling to the left" 
leads to the same final state as "tunneling to the right" . 
We should therefore actually compare AE to 2A£? quar ti c . 



3. Two-level system 

In the present paper we use the energy difference be- 
tween the two lowest eigenvalues of the Schrodinger equa- 
tion to describe coherent quantum oscillations between 
the two classical stationary states of the system. This 
simple approach is only valid as long as the energy bar- 
rier is sufficiently high. In this case the difference between 



the two lowest energy eigenvalues (energy splitting) is 
small. Furthermore, the sum and the difference between 
the eigenfunctions corresponding to the two lowest en- 
ergy eigenvalues are well localized in one of the minima 
of the potential U(Q) and the two lowest energy levels 
arc sufficient to describe coherent quantum oscillations. 

If the barrier becomes too small we need more en- 
ergy levels to describe the quantum dynamics of the 
system sufficiently well. Therefore, we require that 
the difference between the two lowest energy eigen- 
values is significantly smaller than the energy differ- 
ence between the first excited state and the higher en- 
ergy eigenstates. Numerical calculations based on the 
Matematica functions MathieuCharacteristicA and 
MathieuCharacteristicB show that for h w 0.348950 
the difference Ae±2 = £2 — £i between the two lowest 
eigenvalues is the same as the energy difference AS23 = 
£3 — £2 between the first and the second excited state 
and is given by Aei 2 = Ae 2 3 ~ 0.696576. Therefore, we 
have to choose the parameters of our system such that 
h is significantly larger than 0.35. A more detailed anal- 
ysis shows that Ae23 exceeds Aeu by a factor of 2 for 
h > 0.558721, by a factor of 5 for h > 0.936175, and by 
a factor of 10 for h > 1.293984. 

The conditon that h is significantly larger than 
0.35 requires that is significantly larger than 
16. 8 {fvjjp/ EjXj) 2 . As a consequence, we are always in 
a regime where we cannot approximate the system by a 
two-level system if the approximation (|B5|) is applicable. 



Finally we want to mention, that for h f=a 0.329006 
the energy of the first excited state touches the energy 
barrier. For smaller values of h only the lowest energy 
level is below the energy barrier. 
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